Strain-induced pseudo-magnetic field for novel graphene electronics 



Tony Lov^ and F. Guine^^ 

* Network for Computational Nanoelectronics, Hall for Discovery Learning Research 
Purdue University, West Lafayette, IN47907-1791 Indiana, US 
^ Instituto de Ciencia de Materiales de Madrid, CSIC, 
Sor Juana Ines de la Cruz 3, E28049 Madrid, Spain 

Particular strain geometry in graphene could leads to a uniform pseudo-magnetic field of order 
lOT and might open up interesting applications in graphene nano-electronics. Through quantum 
transport calculations of realistic strained graphene flakes of sizes of lOOnm, we examine possible 
means of exploiting this effect for practical electronics and valley tronics devices. First, we found 
that elastic backscattering at rough edges leads to the formation of well defined transport gaps of 
order lOOmeV under moderate maximum strain of 10%. Second, the application of a real magnetic 
field induced a separation, in space and energy, of the states arising from different valleys, leading 
to a way of inducing bulk valley polarization which is insensitive to short range scattering. 

PACS numbers: 



Introduction. The isolation of single layer graphene, 
and the possibility of controlling the electronic carrier 
density T-'T has led to a large research effort, because 
of the novel fundamental features exhibited by graphene 
and also its possible applications [H [5]. Graphene is 
a two dimensional membrane whose electronic proper- 
ties can be controlled by applying a gate voltage. It 
can be shown that the elastic deformations of the mem- 
brane modify the electronic properties, as they play the 
role of an effective gauge field [4, |6l [7], opening a new 
way of tailoring the electronic properties. Strains of 
at least 10% can be induced in graphene without dam- 
aging appreciably its structure [8 . Suspended graphene 
samples show long range deformations on scales of hun- 
dreds of nanometers[9-ll , and strains can be induced in 
graphene samples by different techniques p!2Hl4] . 

Strains can be expected to arise naturally in suspended 
samples [15 . It has been shown that a uniformly varying 
strain leads to a gauge potential which generates an effec- 
tive constant magnetic field [161 Hi] (time reversal sym- 
metry implies that the field has opposite signs in the 
two valleys). The field due to strains can interfere in 
many ways with real magnetic fields [TF. Other schemes 
to use strains [19 to tune the electronic properties have 
also been proposed. Strains can possibly be manipulated 
efficiently in samples with good adhesion to the substrate, 
such as graphene layers grown epitaxially on SiC. Hence, 
the proposed effects in this paper could be tested exper- 
imentally. 

As discussed in[17 , an effective constant magnetic field 
arises from strains varying at a constant rate. We discuss 
here the electronic properties of graphene flakes under the 
combined effect of strains, magnetic fields and disorder. 
We use the setup proposed in [17 as a way of inducing 
an almost constant effective magnetic field. The model 
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FIG. 1: (a) Sketch of an example strain geometry with a max- 
imum strain of 50%. Note that the actual size of the flake used 
in this paper is much larger, i.e. L — W — lOOnm, and with 
much smaller maximum strain (see text). Conductance as 
a function of Fermi energy and strained geometry character- 
ized hy R/W for (b) zigzag and (c) armchair edge ribbons, 
in the absence of real magnetic field and edge disorder. The 
dimension of the graphene flake is L — W — lOOnm. 

and method of calculation are discussed in the next sec- 
tion. Then, we present the main results. The last section 
discusses the main conclusions of our work. 

The model. A strain distribution, Uij^ leads to the 
effective gauge field[6| [7( fTH]. 

where /3 ^ d\og{t)/d\og{a)\a=ao ^ 2 — 3, ao ~ 1.4 A is the 
lattice constant, t ~ 3eV is the nearest neighbor coupling 
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FIG. 2: The non-equilibrium current density sign (j 3,) x |j(a::,y)| 
intensity plot for a 100 x lOOnm graphene flake under a strain of 
R/W = 5 (equivalent to Bs ~ 9T) at ep corresponding to fllling 
factor u = 2 (a) without edge disorder and (b) with edge disorder. 
Similar plots in (c-d) but for fllling factor v = 10 instead. The 
color scale indicated is normalized with respect to max(\}(x,y)\). 
Red/blue color indicates current flowing to the right /left of the 
device, (e) and (f) illustrates the counter-propagating edge states 
a,t V = 2 and = 10 respectively. 



energy and c is a dimensionless constant of order unity. 
The main axes of the graphene lattice in Eq. [l] has been 
chosen so that the x axis coincides with a zigzag direction. 
We assume that the distortions have been induced by 
bending a flake in the way discussed in [17] 
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where R is the bending radius of the deformation ap- 
plied to the flake, see Fig. [l^. This would translates to a 
pseudomagnetic field strength of Bg = cP/aoR = QW/R^ 
where Q is to be determined numerically. 

To compute the transport properties, we adopt a near- 
est neighbor hopping between tt orbitals in the honey- 
comb lattice F20 , 
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where 5dij is the change in bond length, and A is 
the real gauge field, whose effect is incorporated as a 
Peierls phase. We assumed to = 3eV and /3 = 2. 
Semi-infinite leads are assumed for the left and right 
boundaries. The numerical methods that we had em- 
ployed to compute the transmission function T is the 
recursive green function [2 F-'ST and the renormalization 
method [24^. The Landauer formula then gives us the 
zero temperature conductance i.e. 2e^ /h x T[25 . De- 
tails of the implementation had already been documented 
else where [Ml El]. 

We consider square graphene flakes with dimensions 
W = L = lOOnm where the origin for the deformation 
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FIG. 3: Plot of typical energy dispersion as a function of mo- 
mentum along the transport direction for the case of (a) real mag- 
netic fleld B = 9T, (b) pseudomagnetic fleld Bs = 9T and (c) 
Bs = B = 9T. The symbols (g) and O denotes the sign of the 
fleld, while red/black color denotes whether the chiral edge states 
are propagating along the top/bottom edge. The highlighted bulk 
states in (c) are non-chiral i.e. it feels an effective zero magnetic 
fleld. 



in Eq. [2] is the center of the flake. The maximum strain 
exerted would then be along the two edges, which is given 
by mdix{uxx) ^ W/2R. 

Results. We show in Fig{T]3 and c the computed con- 
ductance as a function of R/W and energy for zigzag and 
armchair edge ribbons respectively. Clean quantum Hall 
plateaus are observed for the zigzag ribbons with filling 
factors given by v = 2, 6, 10 . . . = 4n + 2, exactly mimick- 
ing the conventional quantum Hall case. Knowing that 
the first excited Landau energy ei = vrV^HBe in the 
conventional quantum Hall case, we can estimate numer- 
ically the effective magnetic field induced by a deforma- 
tion characterized by the ratio W/R as for Bg = QW/R 
with Q ^ 45T, in good agreement with the estimates 
in [lillT]. 

The effects on the electronic states of a uniaxial strain, 
Uxx{x)^ like the one studied here, depend strongly on the 
orientation of the lattice with respect to the direction of 
the strains. When the y axis coincides with an armchair 
direction, a gauge field along the x direction is generated, 
Ax{y)^ and an effective magnetic field ensues (see for ex- 
ample simpler case of a one-dimensional ripple [28,). On 
the other hand, when the y axis is along a zigzag direc- 
tion, the gauge field due to the uniaxial strain can be 
written as Ay{y). This gauge field does not induce an ef- 
fective magnetic field, and leaves the electronic spectrum 
unchanged. The results in Fig. [l] are in agreement with 
this analysis. Hereafter, we shall only consider zigzag 
ribbons with strain geometry R/W = 5, which yields 
Bg ^ 9T. 



3 




0.0 0.1 0.2 0.3 0.2 0.5 1.0 1.52.0 



ef(eV) A (A) 

FIG. 4: Zero temperature conductance as a function of (a) Fermi 
energy ef and (b) edge disorder A, for different values of A and 
Cf respectively. Zigzag ribbon is assumed and the strain geometry 
employed corresponds to W/R = 5. The solid line for (a) is for the 
case of zero disorder. 



Unlike the case in the Quantum Hah effect, the edge 
states are not protected by time reversal symmetry, and 
they can be affected by elastic backscattering. Fig. |2] a 
and c plots the non-equilibrium current density at Fermi 
energy corresponding to u = 2 and u = 10. For z/ = 2, 
strains induce two edge modes, which propagate in op- 
posite directions. Time reversal symmetry implies that 
these two modes are localized at the same edge, in agree- 
ment with the numerical results. In general, the compres- 
sive strained edge would acquire two more modes than 
the other edge. The zigzag boundary condition used here 
do not mix the K and K' valleys, leading to a clear dis- 
tinction of the edge modes. 

The zigzag boundary condition (where a given edge can 
be characterized by a majority sublattice) is generic for 
graphene edges [29 , except the armchair one. The edge 
modes induced by strains have a characteristic width of 
order of the effective magnetic length associated to the 
strain field, ig ^ \/ {glqL) / {j3u) ^ where L is a typical 
length which describes the variation of the strain, and u 
is the average strain. Typically is ^ <^o, and the contin- 
uum zigzag boundary conditions discussed in [29 describe 
well the numerical results. 

The physical origin of the edge currents is explained 
in Fig. and b, which plots the corresponding energy 
dispersion along the transport direction for the case with 
a real and pseudomagnetic field respectively. For the lat- 
ter, one makes the observation that for a given current 
direction, the edge states on the two edges are valley 
polarized i.e. quantum valley Hall effect. This effect is 
analogous to the quantum spin Hall effect [30 . In both 
cases, the net pseudo-gauge field of the system is zero, 
but finite and opposite for each spin/valley. However, in 
this case, short-range scattering would couples the val- 
leys. Since the counter propagating edge states residing 
along a particular edge belongs to different valleys, inter- 
valley processes lead to backscattering. In the presence 
of edge disorder, substantial backscattering can occurs 
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FIG. 5: (a) shows the conductance as a function of real magnetic 
field and Fermi energy, calculated for non-disordered zigzag ribbon 
with a strain geometry corresponding to W/R = 5, which is equiv- 
alent to Bs ~ 9T. (b) plots the current density at = O-leV for 
the condition B = Bg = 9T, for perfect edge (top) and disorder 
edges (bottom). Similar plots for (c), except now for the condition 
of B = 13T and Bg = 9T. The color scheme employed is the same 
as Fig. [2] 



and Anderson localization spots can be observed (see Fig. 

A more quantitative evaluation of the impact of edge 
roughness follows. 

The edge roughness is generated using a procedure out- 
lined in[27 . The roughness morphologies generated were 
assumed to have an autocorrelation length of 2nm and a 
root-mean-square roughness A, where 0.2A< A < 5 A is 
used as the disorder parameter in this work. Fig. [4^ plots 
the energy resolved conductance for various values of A. 
Disorder at the edges opens a transport gap e^, and the 
flake shows insulating behavior near the neutrality point. 
When the disorder increases, Cg saturates to a value of 
the order of two times the energy of the first effective 
Landau level, modulated by the strain via Cg ex ^/W/R. 
In the nanoribbon counterpart, it was found to behave 
as Cg (X 1/W[31 . The pseudomagnetic length is given by 
is cx y^^R/W^ which approximates the spatial extent of 
the ground state Landau wavefunction from the edges, 
can be of order of lOnm at moderate strain of W/R = 5. 
The estimates in the Supplementary Information in [16] 
suggest a localization length for the edge modes, due to 
intervalley scattering of order ^ ~ (^s/<^o)^^^"^7 where riy 
is the one dimensional density of vacancies at the edges. 
For a rough edge on atomic scales, Uy ~ 10nm~^, so that 
^ ~ W, and localization can be expected for the strains 
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and dimensions used here. 

Next, we examine the interplay between real and pseu- 
domagnetic field. Their combined contribution is such 
that the valleys feel an effective field of A ± A. Fig. 
[5^ shows the conductance for varying B field but at 
Bs ~ 9T. Splitting of the conductance plateaus steps 
from 4n + 2 into 2n + 2 can be observed, with decreasing 
plateau width for n = 0. An interesting situation arises 
when B = Bg, where Fig. [Sj^ illustrates the correspond- 
ing energy dispersion. The cancellation of the field for 
valley leads to a recovery of the bulk band disper- 
sion. Hence, the conductance goes as ~ ce/ + 2, where 
the latter contribution comes from the K valley. This 
provides a feasible avenue for producing valley polarized 
electrons, since it does not depend on good quality at the 
edges like in previous proposal [32 . We note that in gen- 
eral, these bulk states are not protected by time reversal 
symmetry. Fig. [sJd plots the current density at the con- 
dition B = Bg = 9T, with/without edge disorder. Their 
respective conductances (in 2e^ /h) are 5 and 3.3 respec- 
tively. At larger magnetic field, say B = 13T > Bg = 9T, 
edge states emerge and now the conductance remains at 3 
with/without edge disorder. Fig. ^ shows the respective 
current density plots. 

Conclusions. We have analyzed in detail the electronic 
properties of graphene fiakes under the combined effect of 
strains, magnetic fields and disorder. We analyze numer- 
ically the electronic spectrum of fiakes of realistic sizes. 



~ lOOnm with strain distributions that generate a con- 
stant effective magnetic field. Strains induce gaps in the 
bulk spectrum and propagating modes along the edges. 
The edge modes are valley polarized, and strongly sup- 
pressed at armchair edges, where the boundary condition 
strongly mixes the different valleys [33 . The two lowest 
edge modes are localized at the same edge. In general, 
the number of edge modes differs by two between the two 
edges. 

Possible device applications are addressed. First, dis- 
order leads to backscattering, which is more significant 
for the lowest edge (see Fig. ^p)- We find a clear trans- 
port gap near the Dirac energy, leading to new ways of 
developing graphene transistors. Second, the interference 
between real magnetic fields and the gauge field due to 
strains lead to the separation in space and in energy of 
the states from the two valleys. Because of this separa- 
tion, the valley polarization achieved in this way is not 
much affected by intervalley scattering, opening a way of 
obtaining valley polarized bulk currents. 
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